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Kinetic Turbulence 
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Abstract The weak collisionality typical of turbulence in many diffuse astrophysi- 
cal plasmas invalidates an MHD description of the turbulent dynamics, motivating 
the development of a more comprehensive theory of kinetic turbulence. In particu¬ 
lar, a kinetic approach is essential for the investigation of the physical mechanisms 
responsible for the dissipation of astrophysical turbulence and the resulting heating 
of the plasma. This chapter reviews the limitations of MHD turbulence theory and 
explains how kinetic considerations may be incorporated to obtain a kinetic theory 
for astrophysical plasma turbulence. Key questions about the nature of kinetic tur¬ 
bulence that drive current research efforts are identified. A comprehensive model of 
the kinetic turbulent cascade is presented, with a detailed discussion of each compo¬ 
nent of the model and a review of supporting and conflicting theoretical, numerical, 
and observational evidence. 


1 Introduction 


The study of turbulence in astrophysical plasmas has almost exclusively employed 
a magnetohydrodynamic (MHD) description of the turbulent dynamics, treating the 
magnetized plasma as a single fluid, an approximation valid for large-scale, low- 
frequency dynamics in the strongly collisional limit. Yet, the plasmas in a wide va¬ 
riety of turbulent astrophysical environments often violate one or more of the condi¬ 
tions required by the MHD approximation, particularly on the small scales at which 
dissipation mechanisms act to damp the turbulent fluctuations, ultimately leading to 
heating of the plasma. The study of the turbulent dynamics at small scales and of 
the physical mechanisms responsible for the dissipation of the turbulence generally 
requires a kinetic treatment. Thus, it is necessary to leave behind the comfortable 
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surroundings of the theory of MHD turbulence and enter the uncharted territory of 
the evolving theory of kinetic turbulence. 


1.1 Quantitative Characterization of Plasma Turbulence 

Turbulent systems are typically described theoretically by a spectral decomposi¬ 
tion of the broadband spatial fluctuations into a sum of plane wave modes, each 
characterized by its three dimensional wavevector, phase, and amplitude. An en¬ 
ergy spectrum of the turbulent fluctuations therefore provides a useful quantitative 
description of the turbulent system. In a magnetized plasma, the three-dimensional 
wavevector space can be reduced to two dimensions by assuming axial symmetry 
about the direction of the equilibrium magnetic field, requiring only the specification 
of the turbulent power with respect to the cylindrical components of the wavevec¬ 
tor, and k\^y The nature of the dynamics in the different ranges of the kinetic 
turbulent cascade can be quantitatively characterized by two properties: (1) the one¬ 
dimensional magnetic energy spectrum in perpendicular wavenumber, EB{k±)', and 
(2) the wavevector anisotropy, or the distribution of turbulent power in wavevec¬ 
tor space. Here Egik^) is defined such that the total magnetic energy is given by 
Eb = / dk±EB{k±). For Alfvenic turbulence that is driven isotropically at the outer- 
scale wavenumber ko, the conjecture of critical balance implies that the turbulent 
power Alls a region of the cylindrical wavevector space satisfying k|| < 
Specification of the scaling of the boundary of this region, k\\ is sufficient to 

completely characterize the anisotropic distribution of turbulent power. 


1.2 Limits of MHD Treatment of Astrophysical Turbulence 

The limitations of an MHD treatment of astrophysical turbulence can be illumi¬ 
nated by considering the domain of applicability of MHD turbulence theory within 
the broader context of plasma turbulence. Beginning with the general theory of the 
turbulent cascade of kinetic energy in hydrodynamic systems, we consider the modi¬ 
fications required to describe the turbulent energy cascade in the magnetized plasma 
systems relevant to astrophysical environments. 


1.2.1 From Fluid to Kinetic Models of the Turbulent Cascade 

The limitations of MHD turbulence theory can be illustrated most clearly by a qual¬ 
itative comparison of the features of nonlinear cascade of energy in hydrodynamic 
turbulence, MHD turbulence, and kinetic turbulence. 

In hydrodynamic systems, turbulent motions are driven at some large scale L, 
denoted the driving or energy injection scale. Nonlinear interactions serve to trans- 
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fer the turbulent kinetic energy to motions at ever smaller scales, until reaching a 
small scale ly at which dissipation via viscous damping is sufficient to terminate the 
turbulent cascade. For typical hydrodynamic systems, a large dynamic range exists 
between the driving and dissipation scales, L/Zy 3^ 1. In that case, one may define 
an inertial range of scales I within which the effects of the driving and dissipa¬ 
tion are negligible, I ^ ly. Within the inertial range, there exists no particular 
characteristic length scale, so the dynamics of the turbulence in the inertial range 
is found to be self-similar, and a simple application of dimensional analysis is suf¬ 
ficient t o describe accurate ly the steady-state hydrodynamic turbulent cascade of 
energy ( Kolmogorovl 1941 ). A qualitative diagram of the kinetic energy wavenum¬ 
ber spectrum for the hydrodynamic turbulence cascade is shown in Figure[Tla). 

In magnetohydrodynamic systems, the turbulence theory must be modified in 
three important ways. First, the dynamics of two turbulent fields, the velocity and 
the magnetic field, must be described, so the cascade of both kinetic and magnetic 
energy is mediated by nonlinear turbulent interactions. Second, fluctuations of the 
two turbulent fields are dissipated by distinct mechanisms, viscosity for the veloc¬ 
ity and resistivity for the magnetic field. The characteristic length scales of viscous 
dissipation ly and resistive dissipation /,j need not be equal, and their ratio is char¬ 
acterized by the magnetic Prandtl number Pr™ = ly/lri- Third, the magnetic field 
in the plasma establishes a preferred direction, leading to distinct dynamics in the 
direction parallel to the magnetic field and in the plane perpendicular to the mag¬ 
netic field. In addition, the magnetic tension provided by the magnetic field supports 
a type of linear wave, the Alfven wave, which has no counterpart in the hydrody¬ 
namic case, transforming the nature of the turbulent motions from hydrodynamic 
vortices to magnetohydrodynamic waves. This third complication is the most sig¬ 
nificant change from hydrodynamic turbulence, and leads to the inherent anisotropy 
of MHD turbulence, where turbulent energy is transferred more rapidly to small 
perpendicular scales l± than to small parallel scales /||. Nonetheless, despite these 
significant differences, the overall qualitative picture of the turbulent energy cascade 
in MHD turbulence bears a striking resemblance to the hydrodynamic case. 

Consider, in particular, the simplified case of MHD turbulence in a Pr„, = 1 
plasma, so there exists a single dissipation scale I4 = ly = l-q. One may define an 
MHD inertial range, l± ^ l^, directly analogous to the hydrodynamic case. 
Due the anisotropy of the turbulent energy transfer, the turbulent dynamics are op¬ 
timally described with respect to the perpendicular scale l±. The evolution of the 
parallel scale is det ermined in terms of the perp endicular scale by the condition 
of critical balance ( Goldreich and Sridhail 1995h . so that /|| 0 = Z^. The exponent 
q describes the scale-dependent anisotr opy of the MHD turbulent ca scade, where 
q = 213 in the Goldreich-Sridhar model ( Goldreich and Sridhail 1995h . and q=\l2 
in the Boldyrev model dBoldvrevt l2006h . Similar to the hydrodynamic case, in the 
MHD inertial range, there exists no characteristic length scale, so the dynamics of 
MHD turbulence is found to be self-similar as well. Therefore, for the Pr„, = 1 case, 
the MHD turbulence theory appears nearly the same as the hydrodynamic turbu¬ 
lence theory, with a few minor changes: (1) the turbulent cascade is described by 
the perpendicular scale Z^ rather than an isotropic scale Z; (2) there exists a scale- 
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Fig. 1 (a) Wavenumber spectrum for kinetic energy in hydrodynamic turbulence, from the driving 
scale, L, through the inertial range, to the viscous dissipation scale, /y. (b) Perpendicular wavenum¬ 
ber spectrum for total energy E = Ei^+EBva MHD turbulence with Pr,,, = 1, from the driving scale, 
L, through the inertial range, to the viscous and resistive dissipation scale, Ij — ly = l-^. 
















Kinetic Turbulence 


5 


dependent anisotropy due to the parallel scaling /|| 0 = and (3) the exponent p in 
the self-similar power law solution for the one-dimensional energy spectrum E 
may differ quantitatively from the hydrodynamic solution. But the general quali¬ 
tative picture—a self-similar MHD turbulent cascade of energy from the driving 
scale L, through an inertial range, to the dissipative scale —remains essentially 
the same as the hydrodynamic cascade, as is evident by comparing the diagram 
of the wavenumber spectrum for total energy E = £ 1 ^+ Eb in the MHD turbulent 
cascade in Figure [TJb) to the hydrodynamic case in Figure [Tf a). 

In kinetic plasma systems, this simple qualitative model of the turbulent cascade 
changes dramatically due to the existence of three characteristic length scales and 
new physics associated with each of these scales. The three characteristic length 
scales that come into play in typical conditions for turbulent astrophysical plasmas 
are the ion mean free path A;, the ion Larmor radius p,, and the electron Larmor 
radius p^. The MHD approximation requires the following four conditions; 

1. Nonrelativistic conditions, Vts/c <g; 1 

2. Strongly collisional conditions, A,/l ^ 1 

3. Large-scale motions, pi/l <^\ 

4. Low-frequency dynamics, co/£2i<^ 1 

Here = yJlTjnis is the thermal velocitjQ of species i, o is the typical frequency 
of the turbulent fluctuations, and Eli is the ion cyclotron frequency. It is clear that, in 
the MHD approximation, all three of the characteristic scales above are assumed to 
be infinitesimal compared to the typical scale of the turbulent motions, 1. However, 
in astrophysical plasmas of interest, the turbulent dynamics frequently violate con¬ 
ditions (2) and (3) abov^l. Therefore, it is important to examine more closely how 
these characteristic scales enter into the dynamics of the turbulent cascade in astro- 
physical plasmas, leading to a violation of the MHD approximation and requiring 
the transition to a kinetic description of the turbulent dynamics. 


1.2.2 Violation of the MHD Approximation 


Spacecraft measurements of turbulence in the solar wind provide invaluable guid¬ 
ance for the construction of a theoretical model that describes the energy spectrum 
of the kinetic turbulent cascade. Recent measurements of solar wind turbulence with 
unprecedented temporal resolution enable us to probe the turbulent dynamics down 
to the scale of the electron Larmor radius ( Sahraoui_et^ .12009 jKi^ani_et^ , 2009; 


Alexandrova et al. . 20091 : Chen et all 201o[ S^raoui et all 20 Kx Alexandrova et al 

2012h . Therefore. we now have a fairly complete observational picture of the ki¬ 
netic turbulent cascade in the solar wind over a dynamic range of 10^ from the 
the large energy injection scale at L ^ 10® km down to the scale of the electron 


* Here E is expressed in units of energy, absorbing the Boltzmann constant. 

^ Note that condition (4) is not generally independent of condition (3). For MHD Alfven waves, 
the condition co <C may be alternatively written p,//|| ^ where the ion plasma beta is 
j8; = SmiiTi/B^. Thus, if ^(1)> then condition (4) is roughly equivalent to condition (3). 
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Larmor radius at ^ 1 km. From the large body of turbulence measurements in 
the solar wind ( S^mouietahL2009^Kwani_e^^2009^ Alexandrovaet al. , 20091 : 


Chen et all 2010 ; Sahraoui et al. , 2010l:' Alexandrova et al, , 20121) . we can construct 

a general diagram for the perpendicular wavenumber spectrum of the magnetic en¬ 
ergy in turbulent astrophysical plasmas, shown in Figure EJa). It is important to 
emphasize here that, although the general form of the magnetic energy spectrum 
is well established from observations, the interpretation of this spectrum in terms 
of the characteristic plasma scales requires significant input from plasma kinetic 
theory, and many of the features of FigureEJa) remain topics of active research. 

In Figure El a), the plasma turbulence is driven at some large scale L ^ p;. It is 
generally assumed, in the absence of arguments to the contrary, that the turbulence is 
driven isotropically with respect to the magnetic field, so that the perpendicular and 
parallel components of the driving wavevector ko are equal, k||o ^ ko 1 jL. 

If the plasma conditions at the driving scale satisfy the MHD approximation, then 
the large scale end of the turbulent cascade is described by MHD turbulence theory. 
Although the turbulent fluctuations in an MHD plasma may, in general, be com¬ 
posed of a mixture fast, Alfven, and slow waves, observational and numerical evi¬ 
dence suggests that Alfven waves dominate the turbulent dynamics in typical astro- 
physical plasmas (this point is discussed further below). For the Alfvenic turbulent 
cascade, the one-dimensional magnetic energy spectrum as a function of perpen¬ 
dicular wavenumber kx scales a s Fb where the spectra l index is p = —5/3 

in the Goldreich- Sridhar model (IGoldreich and Sridhaii ll99-5 ) or p = —3/2 in the 
Boldyrev model ( Boldvrev 20061) . The wavevector anisotropy of the anisotropic 
Alfvenic cascade scales as k|| where the values for q are given in ^1.2.11 this 
anisotropic cascade of energy through wavevector space is depicted in FigureEtb). 

As the MHD turbulent cascade transfers energy to smaller scales (higher wavenum¬ 
ber), it eventually reaches the one of the characteristic length scales A,, p,, or Pe, at 
which point the MHD approximation is violated. Here we focus on exploring how 
these length scales enter into the model for kinetic turbulence and what effect they 
have on the turbulent dynamics. For typical conditions in astrophysical plasmas, the 
characteristic length scales are ordered by A; > p,- > Pe, so the ion mean free path A; 
is usually reached first. 

The ion mean free path A, characterizes the collisionality for the motion of 
plasma particles parallel to the magnetic flelcil so it must be compared to the par¬ 
allel wavenumber k||. For k|| A,- ^ 1, the plasma is strongly collisional; for k|| A, ^ 1, 
the plasma is weakly collisional. Fluid approximations, such as hydrodynamics or 
MHD, break down for plasma conditions k||A,' > 1, so kinetic theory is formally re¬ 
quired to describe the plasma dynamics in moderately to weakly collisional regimes. 


^ The Lorentz force limits the perpendicular motion of plasma particles to the particle Larmor 
radius. Since typical astrophysical conditions yield p; ^ A,-, the plasma is essentially always colli¬ 
sionless in the perpendicular direction. Note, however, that because plasma particles cannot move 
beyond the Larmor radius in the perpendicular direction from the magnetic field, this embodies the 
large-scale perpendicular motions, /x S> p,-, with a fluid-like behavior, even under weakly colli¬ 
sional conditions. 
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Fig. 2 (a) Perpendicular wavenumber spectrum for magnetic energy in kinetic turbulence, from 
the driving scale, L, through the MHD inertial range to the ion Larmor radius p,, where the turbu¬ 
lent cascade enters the kinetic dissipation range, and down to the electron Lannor radius p^. The 
transition from collisional to collisionless dynamics occurs at (b) Wavevector anisotropy in 
kinetic turbulence, scaling as in the MHD inertial range, in the kinetic dissipation range, 
and k® (no parallel cascade) beyond electron scales. The transition from collisional to collisionless 
dynamics occurs at k||P; ^ 1. 
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As depicted in Figure|2lb), at some point in the MHD inertial range, the parallel 
scales may reach the scale of the ion mean free path, A:|| A, ~ 1, marking the transition 
from collisional dynamics that is well described by MHD at k|| A, ^ 1 to collisionless 
dynamics that requires a kinetic description at k||A, 3> 1. The condition of critical 
balance determines the relation between the parallel and per pendicular wavenum¬ 
bers of strong MHD turbulence ( Goldreich and Sridhai , 1995h . so we may define the 
perpendicular wavenumber kxc that corresponds to the transition of collisionality at 
k||A, ~ 1. This transition in the perpendicular wavenumber spectrum typically oc¬ 
curs at perpendicular scales larger than the ion Larmor radius, k±cPi < 1, as shown 
in Figure 12 a). For perpendicular wavenumbers k± <C kxc, the strongly collisional 
dynamics is well described by MHD, and for k± ^ k±c the weakly collisional dy¬ 
namics require a kinetic description, as depicted in Figure |2 a). 

For the weakly collisional range k± ^ k±c, it has been shown rigorously from 
kinetic theory that the Alfvenic turbulent fluctuations remain essentially fluid in 
nature ( Schekochihin et al. . 20091) . The Alfvenic turbul ent cascade co ntinues to be 
accurately described by the equations of reduced MHD ( StraussL 19761) and remains 
undamped down to the pe rpendicular scale of the ion Larmor radius, k±pi ^ 1 
( Schekochihin et al. . 20091) . Therefore, although the MHD approximation is for¬ 
mally violated at scales k± > k±c, the MHD description of the anisotropic Alfvenic 
cascade remains applicable, regardless of whether the dynamics is collisional or col¬ 
lisionless, for all scales larger than the ion Larmor radius, k±pi <C 1. Therefore, we 
denote the range of scales L ^ /x ^ P; in the kinetic turbulent cascade as the MHD 
inertial range. 


MHD Inertial Range; The range of perpendicular scales from the large scale 
of energy injection to the scale of the ion Larmor radius, l±^ p,-, includ¬ 
ing both collisional and collisionless regimes. 


On the other hand, compressible turbulent fluctuations associated with the MHD 
fast and slow waves in the MHD inertial range require a kinetic description at all 
moderately to weakly collisional scales, k||A,' > 1 or k \ >ku. These m odes are 
damped both collisionally by ion viscosity at k||A, ^ 1 ( Braginskiil 1 965h and col- 
lisionlessly by ion Landau damping at k||A,' :§> 1 ( Barnesl 1966 ). Therefore, it is 
expected that the damped compressible modes will play at most a subdominant role 
relative to the undamped Alfvenic fluctuations in turbulent astrophysical plasmas. 
In the weakly collisional solar wind, for example, compressible fluctuations gen¬ 
erally contribute less than 10% of the turbulent magnetic energy ( Tu and Marschl 
I 995 I: Bruno and Carbona 2005h in the MHD inertial range. Therefore, we turn our 


attention back to the dynamics of the dominant Alfvenic turbulent fluctuations. 

When the Alfvenic turbulent cascade reaches the perpendicular scale of the ion 
Larmor radius, kxP; 1, the MHD description of the Alfvenic fluctuations breaks 
down completely for two reasons. First, finite Larmor radius effects lead to a decou¬ 
pling of the ions from the turbulent electromagnetic fluctuations at perpendicular 
wavenumbers k±pi > 1. The result is that the non-dispersive Alfven wave in the limit 
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k±pj 1 transitions to the dispersive kinetic Alfven wave in the limit k±pi S> 1. The 
dispersive nature of the Alfvenic fluctuations accelerates the rate of the turbulent 
nonlinear energy transfer, leading to a steepening of the magnetic energy spectrum, 
with a break in the spectrum at the k±pi ~ 1, as shown in Figure |2j a). Second, 
collisionless damping of the electromagnetic fluctuations occurs due to the Landau 
resonance with the ions, with a peak in the ion damping rate around k±pi ^ 1. In 
addition, electron Landau damping can also contribute significantly for all scales 
k±pi > 1. The combined effect of the ion and electron collisionless damping can 
lead t o a further steepen ing of the spectrum for scales k±pi > 1 (iHowes et al 


201 it Howes et al. . 201 ih . Finally, the cascade reaches the perpendicular scale of 
the electron Larmor radius, k±Pe 1, where collisionless damping becomes suffi¬ 
ciently strong to terminate the tur bulent cascade, leading to an exponential drop off 


of the magnetic energ y spectrum (ITerrv and et all 12012t lAlexandrova et all 120121 : 


TenBarge et al.Ll2013l) . MHD turbulence theory cannot describe the dispersive wave 


behavior or the dissipation that occurs via kinetic mechanisms at scales k±pi > 1. 
Therefore, we denote the range of scales l± < p, in the kinetic turbulent cascade as 
the kinetic dissipation range. 


Kinetic Dissipation Range: The range of perpendicular scales at or below the 
scale of the ion Larmor radius, l± < p,, where wave dispersion and collision¬ 
less dissipation play important roles. 


1.3 Importance of Kinetic Turbulence 

Of fundamental importance in the study of astrophysical turbulence is to determine 
the pathway by which the energy of turbulent motions is ultimately converted to 
plasma heat. Astrophysical turbulence is generally driven by violent events or insta¬ 
bilities at large scales, but fluctuations are dissipated strongly only at scales of order 
or smaller than the ion Larmor radius. A kinetic turbulent cascade arises to trans¬ 
fer energy via nonlinear couplings from the large energy injection scales, through 
the MHD inertial range, down to the scale of the ion Larmor radius. The turbulent 
fluctuations begin to be damped when the cascade reaches the scale of the ion Lar¬ 
mor radius, marking the entry into the kinetic dissipation range. Since the dynamics 
within the kinetic dissipation range is typically weakly collisional, the dissipation 
of the turbulent electromagnetic fluctuations must be accomplished via collision¬ 
less mechanisms governed by plasma kinetic theory. The energy thus removed from 
the turbulent fluctuations ultimately leads to thermal heating of the protons, elec¬ 
trons, and minority ions in the plasma. The observational signature of astrophysical 
objects depends strongly on the nature of the plasma heating, so to interpret obser¬ 
vational data requires a detailed characterization of the small-scale, kinetic plasma 
turbulence. 
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For example, as matter in an accretion disk spirals slowly into a black hole, it 
converts a tremendous amount of gravitational potential energy into heat. Several 
physic al mechanisms cont r ibute t o this process. First, the m agnetorotational insta¬ 
bility ( Balbus and Hawley , 1991 ; Balbus and Hawle^ 19981) taps free energy from 
the differential rotation of the accretion disk to drive turbulence on the scale height 
of the disk, L ^ H. The turbulence effectively transports angular momentum out¬ 
ward in the disk, enabling accretion disk plasma to fall down the gravitational poten¬ 
tial and mediating the conversion of gravitational potential energy into kinetic and 
magnetic energy of the MHD turbulent fluctuations. The high temperatures charac¬ 
teristic of the plasma in a black hole accretion disk lead to a collisional mean free 
path A, ~ H, so the turbulent dynamics is weakly collisional. A kinetic turbulent 
cascade is responsible for the transfer of turbulent energy through the MHD inertial 
range down to the scale of the ion Larmor radius, where the turbulent electromag¬ 
netic fluctuations are damped via collisionless mechanisms in the kinetic dissipation 
range. An entropy cascade ultimately mediates the final conversion of this turbulent 
free energy into plasma heat. Therefore, the radiation that is emitted from the hot, 
magnetized plasma is a strong function of the black hole properties and of the char¬ 
acter of the small-scale plasma fluctuations, where the plasma heating occurs. To 
interpret observational data from the Chandra X-ray Observatory, for example, one 
must unravel the details of the kinetic turbulent cascade. 

Developing a mature model of the kinetic turbulent cascade is critical to under¬ 
standing the turbulent dynamics of the kinetic dissipation range, the physical mech¬ 
anisms responsible for the damping of the turbulent fluctuations, and the resulting 
heating of the plasma species. The ultimate goal is to develop a predictive capabil¬ 
ity to estimate accurately the heating of the protons, electrons, and minority ions in 
the plasma based on the plasma parameters and the characteristics of the turbulent 
driving. 


2 Key Questions about Kinetic Turbulence 


The unprecedented availability of high temporal resolution solar wind turbulence 
measurements from current spacecraft missions has enabled the observational char¬ 
acterization of the kinetic turbulent cascade from the large scales of energy injection 
down to the scale of the electron Larmor radius. This has spurred the heliospheric 
physics community to engage actively the topic of the turbulence in the dissipa¬ 
tion range of the solar wind, and has engendered considerable controversy about a 
number of significant issues related to the fundamental character of kinetic turbu¬ 
lence. In particular, the nature of both the turbulent fluctuations in this regime and 
the physical mechanisms responsible for their dissipation remains highly contested 
within the scientific community. Four key questions relevant to the study of the dis¬ 
sipation range of solar wind turbulen ce are identified and are discussed at length in 
a forthcoming review (lHowesLl2015[) : 
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1. What are the limits of validity of using a fluid description of the turbulence in the 
dissipation range, and which aspects of the turbulence require a kinetic descrip¬ 
tion? 

2. Are the linear plasma wave properties relevant to the turbulent fluctuations of the 
dissipation range? 

3. What are the characteristic dynamics of the dissipation range fluctuations? 

4. What physical mechanisms are responsible for the dissipation of the turbulent 
fluctuations and the ultimate conversion of their energy to plasma heat? 


Although these significant questions about the nature of kinetic turbul ence re¬ 
main controversial , a promising model of the kinet i c turbulent cascade (IHowes , 
2008t Howes et al. . 20()8at Schekochihin et al. . 20091 : Howes et al. . 2011 ) has been 


developed that appears to be broadly consistent with most observations of solar wind 
turbulence. This model involves an anisotropic cascade of Alfvenic fluctuations be¬ 
ginning as a cascade of Alfven waves in the MHD inertial range and transitioning 
to a cascade of kinetic Alfven waves subject to collisionless damping in the kinetic 
dissipation range. Yet the cascade of energy from large to small scales described 
by this kinetic turbulence model may not explain all of the fluctuations observed in 
the solar wind. For example, fluctu ations can be gene rated by the action of kinetic 
temperature anisotropy instabilities ( Bale et all 2009 ) that are driven by the spheri¬ 
cal expansion of the solar wind, an effect beyond the scope of this model. Plausible 
arguments exist that suggest some of these additional effects may coexist peacefully 
with the kinetic turbulence, proceeding without being significantly affected by or 
significantly affecting the kinetic turbulent cascade. The remainder of this chapter 
aims to describe in detail the model of the kinetic turbulent cascade and to discuss 
the supporting and conflicting theoretical, observational, and numerical evidence. 


3 A Model of the Kinetic Turbulent Cascade 


A basic theoretical model of the kinetic turbulent cascade in astrophysical plasmas 
has been developed with the aim to describe completely the flow of energy from 
the larg e driving scales of the turbulence to its ultimate fate as thermal heat of the 


plasrn a (iHowesl 120081 : iHowes et al.Ll2008aHSchekochihin et al.L 120091: iHowes et al 
201 ih . We present here a brief outline of this model, before delving into a detailed 
description of each component of the model and a discussion of supporting and 
conflicting evidence. 

Violent events or instabilities first drive turbulent fluctuations of the magnetic 
field and plasma at some large scale, generating a mixture of finite amplitude fast, 
Alfven, and slow waves. If the fluctuations are driven isotropically with velocities 
approximately equal to the Alfven velocity in the plasma, a cascade of strong com¬ 
pressible MHD turbulence will mediate the transfer of the turbulent kinetic and 
magnetic energy to smaller scales. The fast waves cascade to smaller scales isotrop¬ 
ically, while the critically balanced Alfven wave cascade produces an anisotropic 
distribution of Alfven and slow wave fluctuations in this collisional part of the MHD 
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inertial range. The parallel scales of the turbulent fluctuations eventually reach the 
ion collisional mean free path, marking the transition from strongly to weakly colli- 
sional dynamics. The compressible fast and slow wave fluctuations suffer collisional 
damping at the moderately collisional scale of the transition, and collisionless damp¬ 
ing at the smaller, weakly collisional scales. The incompressible Alfvenic fluctua¬ 
tions remain undamped through this transition, so the damped fast and slow waves 
are expected to contribute subdominantly to the turbulence compared to the Alfven 
waves. The Alfven waves continue their cascade undamped through the collision¬ 
less remainder of the MHD inertial range until their perpendicular scales reach the 
ion Larmor radius, marking the transition to the kinetic dissipation range. 

The anisotropic Alfvenic fluctuations at this transition transfer energy into a cas¬ 
cade of kinetic Alfven waves at perpendicular scales below the ion Larmor radius. 
In addition, collisionless wave-particle interactions via the Landau resonance with 
the ions lead to a peak in the ion kinetic damping at the ion Larmor radius, dissi¬ 
pating some fraction of the turbulent electromagnetic fluctuation energy. The un¬ 
damped remainder of the turbulent energy continues as a cascade of kinetic Alfven 
waves to smaller perpendicular scales, forming the kinetic dissipation range at all 
scales below the ion Larmor radius. Throughout this range, electron Landau damp¬ 
ing may cause significant collisionless damping of the turbulent fluctuations, with 
the strength of the damping increasing as the perpendicular scale decreases. At the 
perpendicular scale of the electron Larmor radius, the electron Landau damping be¬ 
comes sufficiently strong to terminate the cascade, leading to an exponential decay 
of the turbulent energy spectrum at the electron scale. 

Thermodynamically, the transfer of free energy from the kinetic and magnetic 
energy of the turbulent electromagnetic fluctuations to free energy in velocity space 
structure of the particle distribution functions is not equivalent to irreversible ther¬ 
mal heating of the plasma. Irreversible plasma heating, and the associated increase 
of entropy, ultimately requires collisions. This is accomplished in a weakly colli¬ 
sional plasma by the ion and electron entropy cascades, dual cascades in physical 
and velocity space that drive fluctuations to small enough velocity-space scales that 
arbitrarily weak collisions are sufficient to achieve irreversibility. This final process 
marks the thermodynamic end of the kinetic turbulent cascade, completing the con¬ 
version of large-scale turbulent fluctuation energy to thermal heat of the plasma. 

This model of the kinetic turbulent cascade implies certain answers to the ques¬ 
tions posed in i|2] so we elucidate those answers here; 

1. A fluid description is applicable for all turbulent fluctuations at scales larger than 
the collisional transition, and for the Alfvenic dynamics at all scales larger than 
the ion Larmor radius. The dynamics and kinetic damping of the compressible 
fluctuations at all moderately to weakly collisional scales, and of the Alfvenic 
fluctuations at the scales of the ion Larmor radius and below, require a kinetic 
description. 

2. The properties of the turbulent fluctuations at scales sufficiently below the driving 
scale are related to the characteristics of the linear kinetic plasma waves. 

3. The dissipation range fluctuations are kinetic Alfven waves. 
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4. Ion and electron Landau damping are the physical mechanisms by which the 
turbulent electromagnetic fluctuations are damped, and the ion and electron en¬ 
tropy cascades mediate the irreversible transition of free energy in the particle 
distribution functions to thermal heat. 

In the following sections, we describe in detail all of the facets of this model of 
the kinetic turbulent cascade, providing supporting theoretical, observational, and 
numerical evidence and reviewing findings in conflict with this model. A general 
diagram of the magnetic energy spectrum and the distribution of turbulent power in 
wavevector space is shown in Figured 


3.1 MHD Inertial Range: From Driving Scales to the Collisional 
Transition 


The turbulence in astrophysical environments is typically driven by some external 
mechanism, often a violent event or large-scale instability, that generates plasma 
motions at some large scale, L ^ p,-. This energy injection scale, often denoted the 
outer scale of the turbulence, is an important characteristic of any turbulent astro- 
physical system, and is conveniently parameterized by the wavenumber, ko ^ 1 jL. 
For the investigation of kinetic turbulence, a convenient dimensionless measure of 
the driving scale is the driving wavenumber, koPi, where kop, ^ 1 indicates that the 
turbulence is driven at large scale compared to the ion Larmor radius. It is generally 
assumed, in the absence of arguments to the contrary, that the turbulence is driven 
isotropically with respect to the magnetic field, so that the perpendicular and parallel 
components of the driving wavevector are equal, k||o ^ kj^o ^ ko. 

If the MHD approximation is satisfied for the turbulent dynamics of the plasma at 


the driving scale, then the lar ge scale section of the MHD inertial range is described 
by MHD turbulence theo ry (lSridhM-and_GQ]dreich, 19^; Goldreich and Sridhar 


1995t Gal tier et al. . 2000t Lithwick and Goldreichi 2001 : Boldyreyf 2006h . If the 


amplitude of the driyen turbulent yelocities are comparable to the Alfyen yelocity 
in the magnetized plasma, then a cascade of strong MHD turbulence arises to trans¬ 
fer energy nonlinearly to high er wayenumbers; for smaller amplitudes, weak MHD 
turbu lence will be generated ( Sridhar and Goldreichi 1994 ; Goldreich and Sridhar , 
1995h . Since most turbulent astrophysical enyironments are belieyed to be driyen 


strongly, and weak turbulence eyentually t ransiti ons to strong turbulence as the 
cascade progresses ( Sridhar and Goldreichi 19941) . we focus here on the case of 
strong MHD turbulence. In general, the finite-amplitude turbulent fluctuations may 
be considered to be a mixture of the three propagating MHD waye modes, the in¬ 
compressible Alfyen wayes and compressible fast and slow wayes, as well as the 
non-propagating entropy mode. The nature of the turbulent cascades of these yar- 
ious characteristic fluctuations haye been elucidated by numerical simulations of 
MHD turbulence: the fast wayes cascade isotropically in wayeyector space, while 
the Alfyen wayes, slow wayes and entropy mode fluctuations cascade anisotrop- 
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ically according to the c ondition of critical balance (IMaron and Goldreichl 12001 


Cho and Lazarian , 2003h . 


The fast wave cascade produces an isotropic on e-dimensional magnetic e nergy 
spectrum EB{k) as observed in simulations dCho and Lazarian!l2003b . Two 

competing models exist that describe the n ature of strong MHD turbul ence for 
Alfven waves, th e Goldreich-Srid har model (IGoldreich and Sridhaii Il995h and the 
Boldyrev model (Boldyrev, 2006[) . The magnetic energy spectrum of the Alfvenic 
turbulent cascade is predicted to scale as Ejj k’]^, where the spectr al index is p = 
—5/3 in the Goldreich-Sridhar model ( Goldreich and Sridhai , 1995 ) and p = —3/2 
in the Boldyrev model ( Boldyrev , 2006h . The anisotropy of the Alfvenic cascade, 
for isotropic driving at wavenur nber kn, is given by ku = k\ ‘‘k^^, where q = 2j3 


in the Goldreich- Sridhar model (IGoldreich and Sridhaii 1 19951) . and o' = 1 /2 in the 


Boldyrev model ( Boldyrev . 20061) . The slow waves and entropy modes are pas¬ 
sively cascaded by the Alfven waves, and therefore adopt the same spectrum and 
anisotropic distribution of pow er as the Alfven waves ( Maron and Goldreichl 2001 


Lithwick and Goldreichl 2001 ). For anisotropic turbulent fluctuations with k_L » kn 


the frequencies of the fast wave fluctuations, which scale as O) k, are generally 
much higher than the frequencies of the Alfven and slow wave fluctuations, which 
scale as W kp so the dynamics of the fast wave cascade are expected to decouple 
from the dynamics of the Alfven and slow wave cascades ( Lithwick and Goldreichl 


2001 : Howes et al. . 2012h . 


The turbulent cascade transfers energy nonlinearly to higher wavenumber fluc¬ 
tuations, as dictated by the MHD turbulence theory, until the parallel wavenum¬ 
ber reaches the transition from collisional to collisionless dynamics, k||A,' ^ 1. 
The perpendicular wavenumber, k±c, that corresponds to A:||A, ~ 1, differs for the 
anisotropic Alfven wave cascade and the isotropic fast wave cascade. For the 
anisotropic Alfvenic cascade, the perpendicular wavenumber of this transition is 
given by k±c ko{kQ?ui)^'^^‘<, whereas, for the isotropic fast wave cascade, it is given 
by kj^c k(){k()Xi)^\ or more simply kj^c-Xi 1- Since ^ < 1, this means that the fast 
wave cascade reaches the collisional transition first, at a smaller wavenumber than 
the Alfven wave cascade. 

For many astrophysical plasmas, the transition for both fast and Alfven waves 
occurs within the MHD inertial range, k±cPi < 1- The compressible fast waves, 
slow waves , and e ntropy modes undergo strong collisional damping by ion viscosity 
( Braginskii . 1965h in the moderately collisional conditions at k± ^ k±c- Any energy 
in the compressible turbulent fluctuations that passes through this transition is ex¬ 
pected to be transferred nonline arly to the kinetic counterparts of the MHD fast and 
slow waves ( Klein et^ 2012 h in the we akly collisional conditions at wavenum¬ 
bers k± ^ k±c- ( Schekochihin et al. . 2009h . The Alfven waves are incompressible, 
involving no motions parallel to the magnetic held, so they are essentially unaffected 
by the transition in collisionality, and the Alfven wave cascade continues unabated 
to higher wavenumbers, k± > 
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3.2 MHD Inertial Range: From the Collisional Transition to the 
Ion Larmor Radius 


Critical balance predicts a scale-dependent wavevector anisotropy given by k^/k\^\^ = 
(^x/^o)^ where q < I for either the Goldreich-Sridhar or Boldyrev models. 
Therefore, at perpendicular wavenumbers within the MHD inertial range suffi¬ 
ciently higher than the driving wavenumber, ^ A:o, the Alfvenic fluctuations 
become anisotropic in the sense that kx^ ^ k^^ In the limit of the MHD inertial 
range k±pi 1, the kinetic dynamics of these anisotropic Alfvenic fluct uations is 
described rigorously by the equations of reduced MHD ( StraussL 1976h . and the 
Alfven wave cascade rem ains undamped down to the perpendicular scale of the ion 
Larmor radius, k±Pi ^ 1 ( Schekochihin et al. . 2009ll . It has also been shown that 
the slow wave and Alfven wave casc ades do not exchange energy in the MHD iner¬ 


tial range (ISchekochihin et al.L l2009h . and the fast waves likewise are not expected 


to exchange energy with the Alfven waves due to the mismatch in frequency, as 
discussed in 1 13.11 Therefore, the dynamics of the Alfvenic cascade throughout the 
MHD inertial range is correctly described by the MHD turbulence theory, even at 
the weakly collisional scales, 

^_Lc- 

The magnetic energy spectrum in the solar wind seems to bear this out. Space¬ 
craft measurements in the supe r-Alfvenic sol ar wind are generally interpreted by 
assuming the Taylor hypothesis (ITavlonIl938h . that frequency of measured tempo¬ 
ral fluctuations is directly related to the wavenumber of spatial variations that are 
swept past the spacecraft. At the frequencies / < 0.4 Hz, corresponding to spatial 
scales larger than the ion Larmor radius, the magnet ic energy spectrum in t he solar 
wind has a spectral index of approximately —5/3 dOoldstein et al.L 1995h . appar¬ 
ently consistent with the prediction of the Goldreich-Sridhar theory for strong MHD 
turbulence. It is worth noting, however, that the velo city spectrum was found have a 
spectral index closer to —3/2 dPodesta et al. . 2007h . in conflict with the Goldreich- 
Sridhar model. Recent work on the evolution of the residual energy, Er = Ei^—Eb, in 
MHD turbulence, however, suggests that these spectral indices may indeed be con¬ 
sistent with the Boldyrev theory, and that the difference in the spectral indices of the 
kinetic a nd magnetic energy spectra is an inherent property of the MHD turbulent 
cascade ( Boldvrev et al.L 2011 . 2012 ). 

The cascade of compressible turbulent fluctuations that passes through the col¬ 
lisional transition will suffer modera t e to st rong collisionless damping by the Lan¬ 
dau resonance with the ions dBarnesl Il966h at all higher wavenumbers, k± » kxc- 
The damping of the compressible fluctuations in the moderate to weakly collisional 
regimes at kx > kxc leads to the theoretical prediction that compressible fluctuations 
will play a subdominant role relative to the undamped Alfvenic fl uctuations in tur¬ 


bulent astro physical plasmas. Studies of interstellar scintillation dArmstrong et al 


198ll 1995h show evidence for a power-law spectrum of density fluctuations over 12 
orders of magnitude in the interstellar medium, suggesting that compressible fluc¬ 
tuations are not entirely damped. But it is not possible from remote astrophysical 
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observations to deduce the relative contributions of compressible and incompress¬ 
ible components of the turbulence. 

In situ spacecraft measurements of turbulent fluctuations in the solar wind, how¬ 
ever, allow a direct determination. The entire turbulent cascade in the solar wind, 
including the driving scales, is weakly collisional, A,/L ^ 1, so spacecraft measure¬ 
ments constrain role of compressible fluctuations in collisionless conditions, 
k±c- Measurements show that the turbul ent fluctuations in the MHD inertial range 
appea r to be dominantly incompressible (ITu and MarschLll995tlBruno and Carbone . 
20051) . where the incomp r essible motions have been shown to be Alfvenic in 


nature ( Belcher and Pavia . 1971 ). The compressible fluct uations generally con¬ 


tribute less than 10% of the turbulent magnetic energy (ITu and MarschL 11995 


Bruno and Carbonell2005h . These compressible fluctuations have typically been in¬ 


terpreted as a possible mixture of fast MHD waves and pres sure balanced struc¬ 
tures (BBSs) ( Tu and Marsch . 19^ Bruno and Carbone . 2005 ). wh ere the latter are 
equivalent to non-propagating slow mode fluctuations with kn = 0 (ITu and Marsch , 


1994tlKellogg and Horburvll2005h . N ote, however , that a recent study using a novel 


method of synthetic spacecraft data (iKlein et al.L 120121) suggests that these com¬ 


pressible fluctuations are not associated with the kinetic counterpart of the fast MHD 
wave, but rather consist o f an anisotropic distribution of kinetic slow wave fluctua¬ 
tions (IHowes et al.Ll2012h . Clearly, more investigation of the kinetic physics of com¬ 
pressible turbulent fluctuations in astrophysical environments, including their damp¬ 
ing via collisional and collisionless mechanisms and the resulting plasma heating, is 
needed. Nonetheless, since only a small fraction of the turbulent energy appears to 
be associated with the compressible fluctuations, we focus our attention henceforth 
on the dominant Alfvenic turbulent fluctuations, as they reach the perpendicular 
scale of the ion Larmor radius, k±pi ^ 1. 


3.3 Transition at the Ion Larmor Radius 


Spacecraft measurements of turbulence in the solar wind demonstrate that the —5/3 
scaling of the magnetic energy spectrum in the MHD inertial range breaks at a fre¬ 
quency around / ^ 0.4 Hz, leading to a steeper spectrum at higher frequencies 
in the kinetic dissipation range. Numerous observational studies have attempted 
to correlate the position of the break with a characteristic plasma time or length 
scale, such as the ion cyclotron frequency, the ion Larmor radius, or the ion iner¬ 
tial length ( GpldstdneLak, 19^ : Leamonet al. .ll99'8[ 19991: Leamonet al.L 200] 


Smith et ak . 2001 : Perri et al. . 2010l : Smith et all 20121" Bourouaine et ak . 201 


but contradictory results have been found. Establishing a convincing correlation has 
likely been elusive because three competing effects may contribute to the dynam¬ 
ics at this transition between the MHD inertial range and kinetic dissipation range: 
(1) the transition from non-dispersive to dispersive linear wave physics as the ions 
decouple from the turbulent electromagnetic fluctuations; (2) a peak in the ion ki¬ 
netic damping; and (3) the possible role of kinetic instabilities, such as tempera- 
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ture anisotropy instabilities ( Bale et al. . 2009h . in generating electromagnetic fluc¬ 
tuations at this scale. 

Based on theoretical considerations of the kinetic plasma physics, the kinetic tur¬ 
bulence model presented here predicts that the transition between the relatively well 
understood MHD inertial range and the significantly more controversial kinetic dis¬ 
sipation range occurs at the perpendicular scale of the ion Larmor radius, k \ p, ^ 1 


( Howesl 2008 : Howes et al. . 2008a; Schekochihin et all 2009t Howes et al. . 2011 ). 
The boundary conditions (in wavevector space) for the nonlinear transfer of energy 
into the kinetic dissipation range are given by the nature of the turbulent fluctuations 
at the end of the MHD inertial range. At this transition at k±pi ~ 1, the wavevec¬ 
tor anisotropy of Alfvenic turbulent fluctuations is given by k±/k\\ ^ (koP;)^^ ^ so 
for a s ufficiently large MHD inertia l range, knp, <C 1, this implies k± 3> k|| since 
q < I ( Goldreich and Sridhar , 19951: Boldyrev . 2006h . This significant wavevector 
anisotropy at the transition is suppor ted by multi-spacecra ft measurements of tur¬ 
bulence in the near-earth Solar wind (ISahraoui et alll2010l) . It follows that, beyond 
this transition, the characteristic wavevector of the fluctuations satisfies kxP; > 1 and 
k||p, <c 1; the Alfveni c solution of linear k i netic theory with such a wavevector is the 
kinetic Alfven wave ( Hasegawa and Sat3. 1989 : Stix, 1992 ). Therefore, the Alfven 
waves of the MHD inertial range are predicted to transfer their energy, via nonlinear 
interactions at the transition k \ p ^ 1, t o kinetic Alfven waves ( Learnpn_gL^, 1^98 


Gruzinovi 19981: Leamon et al.l 1999t Ouataert and Gruzinov , 1999t Howes et al 


2008aHSchekochihin et al.[l2009l) . No nlinear gyrok i netic si mulations of this transi¬ 


tion appear to support this hypothesis (iHowes et al.l l2008bh . reproducing the quali- 


tatiye changes in the electric and magn etic field energy sp ectra measured in the solar 


wind at the scale of the spectral break ( Bale et al. . 20051) . 


Another important effect that occurs at the transition at k±pi ~ 1 is a peak in the 


onance with the ions (Leamon et al.. 1998 

Il999: Leamon et all 20001: Howes et al.. 

2008a; Howesl 20081 Schekochihin et al. 

2009i Howes etak. 2011). This ion ki- 


netic damping becomes increasingly strong as the ion plasma beta increases, and 
is generally non-negligible for plasmas with beta of order unity or larger, jS,- > 1, 
leading t o a simifica nt fraction of the dissipated turbulent energy heating the ions 
dHowesl l2oIa Em 1, in approxi mate agreement with empirical estimates of the 
plasma heating in the solar wind ( Cranmer et al. . 2009; Breech et al.l 2009h . Some 
measurements of the magnetic energy spectrum in the dissipation range of the solar 
wind show a significant steepening to a slope of approximate ly —4 at the ion scales . 


20101 ) . 


flattening to —2.8 spectrum further into the dissipation range (ISahraoui et al 
eyidence suggesting significant ion kinetic damping. 

In a steady-state kinetic turbulent cascade, the turbulent energy reaching the tran¬ 
sition at k±pi ~ 1 that is not damped at that scale will carry on, launching a turbulent 
cascade of kinetic Alfyen wayes in the kinetic dissipation range at kxp, > 1. Al¬ 
though the Alfyen and slow waye cascades do not excha nge energy in the MHD in - 


ertial range, they may exchange energy at this transition ( Schekochihin et al. . 20091) . 


so it is possible that the kinetic Alfyen waye cascade can gain energy that is trans¬ 
ferred nonlinearly from compressible fluctuations in the MHD inertial range. 
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3.4 Kinetic Dissipation Range: Between the Ion and Electron 
Larmor Radius 


Although direct spacecraft measurements of the kinetic dissipation range of tur¬ 
bulence in the near-Earth solar wind have been possible for more than a decade, 
the nature of the turbulent fluctuations in this regime remains a controversial topic. 
Characterizing these fluctuations is one of the key goals in heliospheric physics to¬ 
day, especially because the relevant physical dissipation mechanisms that ultimately 
lead to heating of the plasma depend strongly on the nature of the turbulent fluctua¬ 
tions themselves. 

Many early investigations of the dissipation range in solar wind turbulence im¬ 
plicitly assumed that the turbulent fluctuations in the dissipation range are related 
to the linear wave modes in the plasma. Two main hypotheses have been proposed . 


that the tu r bulence is composed of either kinetic Alfven waves ( Leamon_e^ayJ1998 
Gruzinovi 19981 Leamon et al.l 1999t Ouataert and Gruzinovl 1^9t 


2008: 

201 


_ ___ __jHowes et al... 

Schekochihin et al. , 2009 ) or whistler waves (Stawicki et all 200 ll: Gary et al. 


Narita and Garvl 201ol) . Although these two possibilities generally remain 


the leading c andidates, seve r al oth er possibilities have be en suggested: ion Bern¬ 
stein waves ( Sahraoui et al. . 2012h . ion cyclotron waves ( Jian et all 2009ll . non 
propagating pressure balanced structures (BBSs), or inherently nonli near structures. 


partic ularly highly intermittent coherent structures and current sheets (IServidio et al 

201 ih . 


Direct spacecraft measurements of turbulence in the solar wind at the frequen¬ 
cies / > 1 Hz, corresponding to the kinetic dissipation range, provide important 
constraints on the nature of the turbulent fluctuations. A number of recent stud¬ 
ies employing high temporal resolution spacecraft measurements have found a 
nearly power-law scaling of the magnetic energy spectrum between the ion and 
electron scales with a spectral index of approximately —2.8 ( S^iraoui_etaLL 2009; 


Kivani et al. . 20091: Alexandrova et al. . 2009 : Chen et al. . 2010 : Sahraoui et al. . 201ol) 

These observations of the turbulence over the dissipation range scales raise two 
important questions that any model for kinetic turbulence must answer: (1) What 
causes the magnetic energy spectrum to steepen in the dissipation range; and (2) 
Does significant dissipation of the turbulent fluctuations occur between the ion and 
electron scales? 

In the model of the kinetic turbulent cascade, the boundary conditions in wavevec- 
tor space determined by the anisotropic Alfvenic cascade through the MHD inertial 
range suggest that the turbulent energy is transferred nonlinearly to a cascade of 
kinetic Alfven waves in the kinetic dissipation range, as discussed in 113.31 Here we 
describe the properties of the kinetic Alfven wave cascade at perpendicular scales 
below the ion Larmor radius, k±pi > 1. 

Although MHD Alfven waves are non-dispersive, kinetic Alfven waves become 
dispersive due to the averaging of the ion response over the finite ion Larmor radius, 
a physical effect that increasingly decouples th e ions from the electromagnetic fluc- 
tuations with wavevectors satisfying k^pi > 1 (lHollweglll999l : ISchekochihin et al 
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2009h . A useful for mula combining the l i near frequency in the Alfv en and kinetic 
Alfven wave regimes Howes et al, ( 2006 ): Schekochihin et al. ( 2009h is given by 


0) = k\\VA 



( 1 ) 


In addition, the kinetic Alfven wave is significantly compressible, generating a non¬ 
zero parallel magnetic field fluctuation, 5Zj|| , particularly in the limit of low to mod¬ 
erate plasma beta, j3, < 1 ( Hollweg . 19991 : TenBarge et all 2012 ). 

The model for kinetic turbulence predicts the quantitative scaling of the mag¬ 
netic energy spectrum and the wavevector anisotropy for the kinetic Alfven wave 
cascade. The Kolmogorov hypothesis—that the energy transfer rate is constant due 
to local (in wavenumber space) nonlinear interactions—can be used to predict the 
magnetic energy spectrum for the kinetic Alfven wave cascade in the absence of 
significant dissipation. For k±pi ^ 1, the linear wave frequency increases due to 
dispersion, yielding a scaling co k^\k±. This leads to more rapid nonlinear en¬ 
ergy transfer, steepening the magnetic energy spectrum to a predicted scaling Eb < 


-7/3 


when dissipation is neglecte d dBiskmaEetali 1999t QioandLazari^ 2004 


Krishanand^ah^ i 200^IShaikh and Zankl 12005^ iGaltieii l2006t iHowes et ak. 
2008al: Schekochihin et al. . 2009h . Extending the concept of critical balance— 
that the linear wave frequency and nonlinear energy transfer frequency remain 
in balance—to the kinetic Alfven wave regime leads to a predicted wavevec- 

1 /3 I I I I 

tor anisotropy given by k ii k^ (ICho and Lazarianl I2004t IHowes et al.L 12008a : 


Schekochihin et al.L 120091). 


In addition to the effects of wave dispersion, collisionless damping via wave- 
particle interactions can also play an important role in kinetic turbulence for all 
scales k±pi > 1. In addition to the peak in ion Landau damping at k±pi ^ I dis¬ 
cussed in (13.31 electron Landau damping may also play a significant role for all 
scales k\ Pi > 1, becoming in creasingly strong as the perpendicular wavenumber 


increases dHowes et al.Ll200^ . Although early models of the turbulent energy cas¬ 


cade in the kinetic dissipation range suggested that such strong collisionless Lan¬ 
dau damping would lead to an exponential cutoff of the spectrum before reach- 


1 (IHowes et al 


2008al: Podesta et al.. 20101). subsea 

uent solar wind observations called this pre- 

diction into question (Sa 

hraoui et al. 

. 20091: Kivani et al.. 20091: Alexandrova et al.. 

20091 Chen et al.l 2010; 

Sahraoui et al.l 2010l) and recent kinetic numerical Simula- 


tions have demonstrated that this idea is incorrect (IHowes et al.Ll20111) . 


In addition to collisionless damping via the Landau resonance, if the kinetic 
Alfven wave frequency reaches the ion cyclotron frequency, co Qi, collision¬ 
less damping may occur via the cyclotron resonance with the ions. However, the 
very large MHD inertial range typical of astrophysical plasma turbulence leads to 
highly anisotropic fluctuations at small scales, kj| <C k±, so the kinetic Alfven wave 
frequency typically remains very small compared to the ion cyclotron frequency, 
CO Qj. Therefore, ion cyclotron damping is not predicted to play a strong role in 
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the dissipation of astrophysical turb ulence (Howes et al.Ll2008ah . with a few excep¬ 
tions, such as the inner heliosphere (lHowesll201 It) . 

There is significant evidence accumulating in support of a kinetic Alfven wave 
cascade at the perpendicular scales between the electron and ion Larmor radius, 
but there also remains observational evidence that appears to be unexplained by 
this model. The scaling predictions for the kinetic Alfven wave cascade in the ab¬ 
sence of dissipation have been corroborated by simulations using electron MHD, 
a fluid limit which describes the dynamics of kinetic Alfven waves in the limit 
k\\ <C k±, but does not resolve the physics of collisionless dissipation. Specifically, 


these simulations reproduce the predicted magnetic energy scaling, Eb 




( Biskmap et ah . 1999t Cho and Lazarian , 2004 , 2009 : Shaikh and Zank . 2009), and 


wavevector anisotropy, kn ( Cho and Lazarian , 2004 2009). The magnetic en¬ 


ergy spectrum from these fluid simulations, ho wever, is not consistent with the oh 
served spectral index of approximately —2.8 ( Sahraour^etakL 2009Kiyani et al. 


2009t Alexandroya et al.. 2009t Chen et al.. 2010HSahraoui et ^ 2010l). 

It has been recently suggested that the combined effects of collisionless dissipa¬ 
tion and nonlocal energy transfer can lead to a further steepening of the magnetic 


energy spectrum beyond k , for scales k±pi > 1 (iHowes et al.L 1201 Ih . For a hy 


drogenic plasma of protons and electrons, the dynamic range between the ion and 
electron Larmor radius for unity temperature ratio is p,/pe ~ 43, so there is little 
room for an asymptotic range of perpendicular scales satisfying the requirements 
1/p; <C ^ 1 /pe- Therefore, it should come as no surprise that a self-similar 
spectrum with a spectral index of —7/3 is not obseryed—throughout the range of 
perpendicular scales between the ion and electron Larmor radius, the transition at 
kxPi ^ 1 and strong kinetic dissipation at k^Pe ^ 1 may significantly affect the 
turbulent dynamics. 

Nonlinear gyrokinetic simulations of turbulence in the kinetic dissipation range 
seem to bear this out. A simulation oyer the entire range of scales from the ion 
to the electron Larmor radius, which resolyes the physics of collisionless ion and 
electron damping, produces a nearly power-law magnetic energy spectrum with a 
spectral index of —2.8, in remark able quantitatiye agreement with the solar wind 
measurements ( Howes et al. . 201 ih . Additional gyrokinetic simulations support the 


1/3 A 

predicted scaling of t he wayeyector anisotropy, k|| k^ (ITenBarge and Howes . 
2012 : TenBarge et all 2013 ). 

Direct spacecraft measurements of dissipation range turbulence in the solar wind 
haye yielded other lines of eyidence in support of or in conflict with the model of a 
kinetic Alfyen waye cascade. A k-flltering analysis of multi-spacecraft data from the 
Cluster mission establishes that the plasma-frame fluctuation frequencies are consis¬ 
tent with linear dispersio n relation of the kinetic Alfyen waye and inconsistent with 
that of the whistler waye ( Sahraoui et al. . 2010l) . A study combining measurements 
of the ratio of electric to magnetic held fluctuation amplitudes and of the magnetic 
compressibility haye shown that the small-scale fluctuations agree well with pre- 
dictions for kinetic Alfyen wayes and are inconsistent with that for whistler wayes 
( Salem et all 2012 ). An examination the compressibility of turbulent fluctuations in 
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the weakly collisional plasma in the MHD inertial range finds evidence of negligi¬ 
ble energy in the fast wave mode, suggesting that all large-scale turbulent energy is 
transferred, via the anisotropic Alfv enic cascade, to kinetic Alfven waves, w ith little 
energy coupling to whistler waves ( Howes et all 2012t Klein et al. . 2012 ). Investi¬ 
gations of the magnetic helicity of turbulent fluctuations as a function of the angle 
of the wavevector with respect to the loca l magnetic field direction finds a broad re - 
gion of positive helicity at oblique angles ( He_et^ ^ 20nt Podesto and Garvi 2011 ), 
as expected for kinetic Alfven waves ( Howes and Ouataerll 201ol) . but a small re¬ 
gion corresponding to nearly parall el wavevectors that is consistent with eit her ion 


cyclotron waves or whistler waves (IHe et al.Ll201 ItiPodesta and Garvil201 Ih . 


The presence of either ion cyclotron or whistler waves with nearly parallel 
wavevectors is not explained by the model for kinetic Alfven wave turbulence, but 
these fluctuations may be driven by the action of kinetic temperatu re anisotropy 


instabilities in the spherically expanding solar wind ((Bale et al 


20091) . These i 


insta¬ 
bilities typically generate relatively isotropic fluctuations (with respect to the local 
mean magnetic field direction), having wavevector components ^ k||P, ^ 1. 
Since the anisotropic Alfvenic cascade produces fluctuations with k|| ^ k±, and 
since Alfvenic frequencies in equation ([T]i scale linearly with the parallel com¬ 
ponent, 0 ) k||, these anisotropy-driven fluctuations are expected to have a much 

higher frequency and to occupy a different regime of wavevector space than turbu¬ 
lent fluctuations of the turbulent cascade. Therefore, it is reasonable to expect that 
any kinetic-instability-driven fluctuations may persist without significantly affect¬ 
ing, or being significantly affected by, the turbulent fluctuations of the anisotropic 
Alfvenic cascade. 

In conclusion, although the nature of the kinetic turbulence at the perpendic¬ 
ular scales between the electron and ion Larmor radius has not been established 
conclusively, there appears to be significant evidence for an anisotropic cascade of 
kinetic Alfven waves in the kinetic dissipation range. Collisionless dissipation via 
the Landau resonance with the electrons appears to be play a non-negligible role 
in steepening the magnetic energy spectrum beyond the dissipationless prediction. 
But this damping is not strong enough to halt the cascade, so the kinetic turbulence 
continues down to the perpendicular scale of the electron Larmor radius, at which 
point strong kinetic dissipation can effectively terminate the turbulent cascade. 


3.5 Kinetic Dissipation Range: Termination at Electron Larmor 
Radius 

Ultimately, the kinetic turbulent cascade reaches the perpendicular scale of the elec¬ 
tron Larmor radius, k±Pe ^ 1. At this scale, the linear collisionless damping rate due 
to electron Landau damping reaches a value y/o) ^ 1, sufficiently strong that the tur¬ 
bulent magnetic energy cascade is terminated. A simplified analytical treatment of 
the turbulent cascade undergoing this dissipation suggests the spectrum will develop 
an exponential fall-off Eg kj^’^exp(—kxPe) setting in at the perpendicular scale 
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of the electron Larmor radius, k^pe 1 (ITerrv and et all 12012h . As the amplitudes 


of the turbulent fluctuations are diminished by damping, the strong kinetic Alfven 
wave turbulence eventually drops belo w critical balance a nd becomes weak dissi¬ 
pating kinetic Alfven wave turbulence ( Howes et al. . 201 ih . It has been conjectured 
that the transition back to weak turbulence leads to an inhibition o f the parallel 


cascade, so the parallel number of the fluctuations remains constant (iHowes et al 


201 ih . as shown in in Figure|2b). 


A recent study of a sample of 100 solar wind magnetic energy spectra at the 
electron scales shows that all of these spectra may be fit by an empiri cal form 


Eb k“exp(—kxPe), where —2.5 > Ct > —2.8 (lAlexandrova et all 120121) . A non¬ 


linear gyrokinetic simulation of the turbulence over the range 0.12 < k^Pe < 2.5 
yields an energy spectrum demonstrating an exponential fall-off that is quantita¬ 
tively fit by Fg oc kj^'^expf—k I Of,), fu rther supporting the model of the kinetic 


turbulent cascade (TenBarge et akf 20131) . A refined model of the turbulent cascade 


(iHowes et al.L 120111) —incorporating the weakening of the nonlinear turbulent en¬ 


ergy transfer due to dissipation, the effect of nonlocal energy transfer, and the linear 
collisionless damping via the Landau resonance—fits the shape of the spectrum 
well. This provides compelling evidence that collisionless damping is the dominant 
mechanism for the dissipatio n of the kinetic turbul ent cascade, marking the end of 


the kinetic dissipation range (ITenBarge et al.Ll2013h . 


A number of recent works have suggested instead that dissipation in current 
sheets is the dominant dis sipation mechanism for plasma turbulence, based on fluid 
simulations usin g MHD dOmitruk et al. . 20041: Servidio et al. . 2009 . 2010l 2011 ) 
and Hall MHD ( Dmitruk and Mat thaeusi l2006l). hybrid simulations with kinetic 
ions and fluid electro ns ( Parashar et al.L 2009 , 2010l: Markovskii and Vasauez , 2011 


Servidio et al.L 2012h. and observatio nal studies of large-scale discontinuities in the 


solar wind ( Osman et ah . 201 1[ 2012h . However, all the numerical work upon which 
this conclusion has been based employ a fluid description of the electrons which 
does not resolve the dominant collisionless dissipation mechanism of Landau damp¬ 
ing. In a collisionless plasma, current sheets supporting small-scale reconnection 
with a guiding magnetic field are expected to form at the perpendicular scale of the 
electron Larmor radius, k±Pe ~ 1 ( Birn and PriestLEoObh . A recent gyrokinetic sim¬ 
ulation over the range of electron scales 0.12 < k±pg < 2.5, which resolves both col¬ 
lisionless electron Landau damping and the formation of current sheets at k±Pe ^ 1, 
finds dissipati on via current sheets to be su b-dominant compared to linear collision¬ 
less damping ( TenBarge and Howes . 20131) . 

These results establish fairly secure observational and numerical grounds that 
the electromagnetic fluctuations of the kinetic turbulent cascade are dissipated at 
the perpendicular scale of the electron Larmor radius, k^Pe ~ 1. Further work is on¬ 
going to identify the dominant physical mechanisms for the dissipation of these tur¬ 
bulent electromagnetic fluctuations. Although this dissipation terminates the kinetic 
dissipation range of electromagnetic fluctuations at electron scales, there remains 
the final matter of identifying the physical mechanism mediating the conversion 
of this turbulent fluctuation energy irreversibly to thermodynamic heat in a weakly 
collisional astrophysical plasma. 
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3.6 Irreversible Heating Via the Ion and Electron Entropy 
Cascades 


At the perpendicular scales of the ion and electron Larmor radius, collisionless 
wave-particle interactions via the Landau resonance damp the turbulent electromag¬ 
netic fluctuations. In the absence of collisions, this process conserves a general¬ 
ized energy, generating nonthermal struct ure in velocity space of the correspond¬ 
ing plasma particl e distribution functions ( Howes , 2008 ; Schekochihin et all 20091 : 
Plunk et ^ 2010l) . Boltzmann’s //-theorem dictates that, in a kinetic plasma, col¬ 
lisions are required t o increase the entro py and therefore achieve irreversible ther¬ 
modynamic heating ( Howes et al. . 20061) . In the weakly collisional plasmas of as- 
trophysical environmen t s, an entropy cascade —a nonlinear phase mixing process 
(IDorland and Hammetd Il993h that drives a dual cascade in physical and veloc¬ 
ity space—mediates the transfer, at sub-Larmor radius scales, of the nonthermal 
free energy in the particle distribution functions to sufficiently small scales in ve¬ 
locity space that arbitrarily weak collisions can manifes t irreversibility, increasing 


the entropy and th ermodynamically heating the plasma ( Schekochihin et al. . 20091 : 


Plunk et alll2010l) . This inherently kinetic physical mechanism represents the final 


element of the kinetic turbulent cascade, governing the final transition of the turbu¬ 
lent energy to its ultimate fate as plasma heat. 

The ion entropy cascade in two-dimensional plasma systems (in the plane per- 


reticallv (Schekochihin et all 20091 

Plunk et all 2010b and verified in gvrokinetic 

numerical simulations (ITatsuno et al. 

. 2009; Plunk and Tatsuno 

201 ll). In the inher- 

ently three-dimensional system of Alfvenic plasma turbulence ( 
the effects of the ion entropv cascade in generating structure 
comnonent of velocitv snace (Howes, 20081) and in manifesting 
ical scales well below the peak in the collisionless ion damping 
have been identified numerically. Yet, a thorough analysis of the 

Howes et all 2011). 

in the perpendicular 
ion heating at phvs- 

(Howes et all 201 lb 

ion and electron en- 


tropy cascades in kinetic turbulence remains to be undertaken. 


4 Conclusion 

Turbulence is found ubiquitously throughout the universe, playing a governing role 
in the conversion of the energy of large-scale motions to astrophysical plasma heat. 
Extending our understanding of astrophysical turbulence from the limited theory 
of MHD turbulence to the more comprehensive theory of kinetic turbulence opens 
up the possibility of ultimately achieving a predictive capability to determine the 
plasma heating due to the dissipation of turbulence. This chapter has outlined a 
theoretical model of the kinetic turbulent cascade describing the flow of energy from 
the large driving scales, through the MHD inertial range, to the transition at the ion 
Larmor radius, and into the kinetic dissipation range, where the energy is ultimately 
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converted to plasma heat. Although significant progress has already been made, 
much research remains to be done to refine the kinetic turbulence model and test its 
predictions using numerical simulations and observational data. 

Since kinetic turbulence includes a number of the physical processes that are in¬ 
herently kinetic—such as collisionless wave-particle interactions and the entropy 
cascade—kinetic numerical simulations will play an essential role in testing the 
predictions of this model. The higher dimensionality of kinetic systems—in gen¬ 
eral, requiring three dimensions in physical space and three dimensions in velocity 
space—demands a huge investment of computational resources to perform numer¬ 
ical simulations. It is tempting to reduce the dimensionality in physical space to 
two-dimensions to lower the computational costs, but doing so fundamentally lim¬ 
its the applicability of the results to turbulent astrophysical plasmas. The reason is 
because the anisotropic Alfvenic turbulence dominating the kinetic turbulent cas¬ 
cade is inherently three-dimensional in physical space: the dominant nonlinearity 
responsible for the turbulent cascade requires both dimensions perpendicular to the 
magnetic field, and Alfvenic fluctuations require variation in the parallel dimension 


(iHowes et al.Ll201 Ih . Therefore, kinetic simulation results can only be directly com¬ 


pared to astrophysical systems if physical space is modeled in three dimensions. 

Observational tests of the kinetic turbulence model should exploit intuition 
from kinetic plasma theory to unravel the dependence of the turbulent proper¬ 
ties on the plasma parameters. The suitably normalized MHD linear dispersion 
relation depends only on two dimensionless parameters, ®iy[iq£) = ®/(kvA) = 
a)pylHD(i3,0); the plasma beta i3, an d the angle between the wavevector and the 
magnetic field 9 ( Klein et ^ 2012 ). In contrast, the linear physics of Vlasov- 
Maxwell kinetic theory depends on five dimensionless parameters, ft>vM = ®= 
0)\/j^{k\\Pi,k±Pi, Pi, TijTe^Vtijc): the normalized parallel wavenumber kj|P,, the nor¬ 
malized perpendicular wavenumber k±pi, the ion plasma beta )3,, the ion-to-electron 
temperatur e ratio Tj/Tj,, and the ra t io of the ion thermal velocity to the speed of 
light Vti/c ( Stix . 1992 : Howes et al. . 2006h ^The dynamical behavior of the kinetic 
plasma—for example, the frequencies, collisionless damping rates, and eigenfunc¬ 
tions of the fluctuations—varies as these dimensionless parameters are changed. 
Therefore, since the kinetic plasma physics depends on these parameters, observa¬ 
tional investigations of kinetic turbulence should strive to analyze measurements 
in terms of the ion plasma beta jS,-, the ion-to-electron temperature ratio Ti/Tg, and 
length scales normalized to a characteristic plasma kinetic length scale, such as the 
ion Larmor radius or ion inertial length. 

Finally, we conclude this chapter on kinetic turbulence with a schematic diagram 
that depicts the salient features of the kinetic turbulent cascade, as shown in Fig¬ 
ure 0 for the case of turbulence in the near-Earth solar wind (the relevant panel is 
noted in parentheses in the description below). The turbulence is driven isotropi¬ 
cally at a large scale, corresponding to a normalized wavenumber koP; 10^^ (a,b). 


In the limit that the turbu lent astrophysical fluctua t ions satisfy the gyrokinetic approxima¬ 
tion, k\\ and a iFrieman and Cheil 1 19821 : iHowes et all l200d : ISchekochihin et all 

l2009l) . the linear phys ics depends on only three dimensionless parameters, = (o/{k\p!A) = 
owes et 'in i2006h . 
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Fig. 3 (a) Diagram of the magnetic energy spectrum and ion entropy cascade in kinetic turbulence, 
(b) Anisotropic distribution of power in (A:^,A:||) wavevector space in kinetic turbulence. 
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Nonlinear interactions serve to transfer the energy from this low driving wavenum¬ 
ber to higher wavenumber through the MHD inertial range, generating a magnetic 
energy spectrum scaling as Eg <== where p = —5/3 or —3/2 (a). This cascade 
is anisotropic in wavevector space, such that turbulent fluctuations fill the shaded 
region below where ^ = 2/3 or 1/2 (b). The Alfvenic turbulence transi¬ 

tions from the MHD inertial range to the kinetic dissipation range at a perpendicular 
wavenumber k±pi ^ 1, a scale at which collisionless ion Landau damping peaks (a). 
The energy transferred via wave-particle interactions to the ion distribution function 
feeds the ion entropy cascade at wavenumbers k±pi > 1, a dual cascade in physical 
and velocity space that mediates the transfer of nonthermal structure to sufficiently 
small scales in velocity space that weak collisions can thermalize the energy (a). The 
remaining turbulent energy that is not collisionlessly damped at ~ 1 is trans¬ 
ferred nonlinearly to a kinetic Alfven wave cascade in the kinetic dissipation range, 
kxP; > 1, leading to a magnetic energy spectrum Eb (a) and a wavevector 

anisotropy k\\ k^ (b). In addition, electron Landau damping becomes stronger as 
the wavenumber increases over the entire range kxp, > 1 (a). Finally, electron Lan¬ 
dau damping becomes sufficiently strong to terminate the kinetic turbulent cascade, 
leading to an exponential fall-off of the magnetic energy spectrum at kxPe ~ 1 (a), 
inhibiting the transfer of energy to higher parallel wavenumber (b), and possibly 
launching an electron entropy cascade (not shown). 
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